Suppression of High-px Neutral Pion Production in Central Pb+Pb Collisions at 
v /snn = 17.3 GeV Relative to p+C and p+Pb Collisions 

M.M. Aggarwal, 1 Z. Ahammed, 2 A.L.S. Angelis, 3 '* V. Antonenko, 4 V. Arefiev, 5 V. Astakhov, 5 V. Avdeitchikov, 5 
T.C. Awes, 6 P.V.K.S. Baba, 7 S.K. Badyal, 7 S. Bathe, 8 B. Batiounia, 5 C. Baumann, 8 T. Bernier, 9 K.B. Bhalla, 10 

V.S. Bhatia, 1 C. Blume, 8 D. Buchcr, 8 H. Biisching, 8 L. Carlen, 11 S. Chattopadhyay, 2 M.P. Decowski, 12 
H. Delagrange, 9 P. Donni, 3 M.R. Dutta Majumdar, 2 K. El Chenawi, 11 A.K. Dubey, 13 K. Enosawa, 14 S. Fokin, 4 
V. Frolov, 5 M.S. Ganti, 2 S. Garpman, 11 '* O. Gavrishchuk, 5 F.J.M. Geurts, 15 T.K. Ghosh, 16 R. Glasow, 8 

B. Guskov, 5 H. A.Gustafsson, 11 H. H.Gutbrod, 17 I. Hrivnacova, 18 M. Ippolitov, 4 H. Kalechofsky, 3 
R. Kamermans, 15 K. Karadjev, 4 K. Karpio, 19 B. W. Kolb, 17 I. Kosarev, 5 '* I. Koutcheryaev, 4 A. Kugler, 18 
P. Kulinich, 12 M. Kurata, 14 A. Lebedev, 4 H. Lohner, 16 L. Luquin, 9 D.P. Mahapatra, 13 V. Manko, 4 
M. Martin, 3 G. Martinez, 9 A. Maximov, 5 Y. Miake, 14 G.C. Mishra, 13 B. Mohanty, 2 * 13 M.-J. Mora, 9 
D. Morrison, 20 T. Mukhanova, 4 D. S. Mukhopadhyay, 2 H. Naef, 3 B. K. Nandi, 13 S. K. Nayak, 7 T. K. Nayak, 2 

A. Nianine, 4 V. Nikitine, 5 S. Nikolaev, 5 P. Nilsson, 11 S. Nishimura, 14 P. Nomokonov, 5 J. Nystrand, 11 
A. Oskarsson, 11 I. Otterlund, 11 S. Pavliouk, 5 T. Peitzmann, 15 D. Peressounko 4 V. Petracek, 18 S.C. Phatak, 13 
W. Pinganaud, 9 F. Plasil, 6 M.L. Purschke, 17 J. Rak, 18 M. Rammler, 8 R. Raniwala, 10 S. Raniwala, 10 N.K. Rao, 7 
F. Retiere, 9 K. Reygers, 16 G. Roland, 12 L. Rosselet, 3 I. Roufanov, 5 C. Roy, 9 J.M. Rubio, 3 S.S. Sambyal, 7 
R. Santo, 8 S. Sato, 14 H. Schlagheck, 8 H.-R. Schmidt, 17 Y. Schutz, 9 G. Shabratova, 5 T.H. Shah, 7 I. Sibiriak, 4 
T. Siemiarczuk, 19 D. Silvermyr, 11 B.C. Sinha, 2 N. Slavine, 5 K. Soderstrom, 11 G. Sood, 1 S.P. S0rensen, 20 
P. Stankus, 6 G. Stefanek, 19 P. Steinberg, 12 E. Stenlund, 11 M. Sumbera, 18 T. Svensson, 11 A. Tsvctkov, 4 
L. Tykarski, 19 E.C.v.d. Pijll, 15 N.v. Eijndhoven, 15 G.J. v. Nieuwenhuizen, 12 A. Vinogradov, 4 Y.P. Viyogi, 2 
A. Vodopianov, 5 S. Voros, 3 B. Wyslouch, 12 G.R. Young 61 ' 2 < 3 < 4 - 5 - 6 < 7 - 8 - 9 < 10 < n > 12 - 13 ' 14 ' 15 ' 16 ' 17 > 18 < 19 > 20 > 21 - 22 

1 (WA98 Collaboration) 
21 University of Panjab, Chandigarh 160014, India 
32 Variable Energy Cyclotron Centre, Calcutta 700064, India 
J,i University of Geneva, CH-1211 Geneva 4, Switzerland 
54 RRC "Kurchatov Institute", RU- 123 182 Moscow 
65 Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 
76 Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6372, USA 
87 University of Jamrau, Jammu 180001, India 
98 University of Miinster, D-48149 Munster, Germany 
109 SUBATECH, Ecole des Mines, Nantes, France 
1110 University of Rajasthan, Jaipur 302004, Rajasthan, India 
1211 University of Lund, SE-221 00 Lund, Sweden 
1312 mT Cambridgej M a 02139 

J ' i13 Institute of Physics, Bhubaneswar 751005, India 
1514 University of Tsukuba, Ibaraki 305, Japan 
1615 Universiteit Utrecht/ NIKHEF, NL-3508 TA Utrecht, The Netherlands 
1716 KVI, University of Groningen, NL-9747 AA Groningen, The Netherlands 
1817 Gesellschaft fur Schwerionenforschung ( GSI), D-64220 Darmstadt, Germany 
191S Nuclear Physics Institute, CZ-250 68 Rez, Czech Rep. 
2019 Institute for Nuclear Studies, 00-681 Warsaw, Poland 
2120 University of Tennessee, Knoxville, Tennessee 37966, USA 

22 * Deceased 
(Dated: June 17, 2008) 

Neutral pion transverse momentum spectra were measured in p+C and p+Pb collisions at y / %N = 
17.4 GeV at mid-rapidity (2.3 < ryiab < 3.0) over the range 0.7 < pt ^ 3.5 GeV/c. The spectra are 
compared to n° spectra measured in Pb+Pb collisions at ^/snn = 17.3 GeV in the same experiment. 
For a wide range of Pb+Pb centralities (N par t ^ 300) the yield of 7r°'s with pt ^ 2 GeV/c is larger 
than or consistent with the p+C or p+Pb yields scaled with the number of nucleon-nucleon collisions 
(JVcoii), while for central Pb+Pb collisions with jV pa rt ^ 350 the ty° yield is suppressed. 

PACS numbers: 25.75.Dw 



The study of hadron production at high transverse and in particular, to detect the possible formation of a 
momentum (pr) is a sensitive tool to characterize the quark-gluon plasma (QGP), i.e., a thermalizcd phase in 
matter created in ultrarelativistic heavy-ion collisions, which quarks and gluons are the relevant degrees of free- 
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dom. Particles at high px result from quark and gluon 
scatterings with high momentum transfer ( "hard scatter- 
ing") which can be described by perturbative quantum- 
chromodynamics (pQCD). The scattered quarks and glu- 
ons will traverse the created medium and fragment into 
the observable hadrons. High-px particle production in 
nucleus-nucleus (A+A) collisions was predicted to be 
suppressed pj as a consequence of the energy loss of 
the scattered partons in the dense matter ("jet quench- 
ing"). Such suppression was observed by experiments 
at the Relativistic Heavy Ion Collider (RHIC) in central 
Au+Au and Cu+Cu collisions at a center-of-mass en- 
ergy of up to v /IiS? = 200 GeV @, i |, H, @. Within 
jet-quenching models the suppression can be related to 
medium properties, such as the initial gluon density or 
the transport coefficient @, S, @| • 

A crucial test for the idea of parton energy loss in the 
hot and dense medium created in A+A collisions is the 
measurement of the ^/snn dependence of high-px hadron 
production @, H, 0, El, El, El . In central Pb+Pb col- 
lisions at the CERN SPS energy of y/s^ = 17.3 GeV 
(piab = 158 ^4GeV/c) the initial energy de nsity , as esti- 
mated from the measured transverse energy 14| , is above 
the critical value e c ~ 1 GeV/fm 3 [TH for the transition 
to the QGP. On the other hand, the initial gluon den- 
sity and the lifetime of a deconfmed phase produced at 
•y/sNN = 17-3 GeV will be significantly reduced as com- 
pared to RHIC energies. Results at SPS energies thereby 
provide a sensitive test of jet quenching model predic- 
tions. 

Results on high-px particle production in central 
Pb+Pb collisions at the CERN SPS have already been 
published 1|| 17 1. However, the interpretation of these 
data has been complicated by the lack of reference p+p 
data to allow to quantify nuclear effects. Instead, param- 
eterizations of p+p data have been employed to search 
for nuclear effects but with substantial system- 

atic uncertainties. Moreover, hadron suppression due 
to parton energy loss might be compensated by an en- 
hancement due to multiple soft scatterings of the incom- 
ing partons prior to the hard scattering process ("nu- 
clear fcx-enhancement" or "Cronin effect"). Measure- 
ments in p+A and d+A collisions at different ^/snn sug- 
gest that such enhancement is significantly stronger at 
v /5ii" m 20 GeV than at = 200 GeV [U]. 

In this letter ir° spectra are presented from p+C and 
p+Pb collisions at ^snn = 17.4 GeV (pi a b = 160 GeV/c) 
measured in the WA98 experiment. Since the nuclear 
fcx-enhancement in p+C is expected to be small [l9|] 
this measurement should provide a useful substitute for 
a p+p reference. Information on the magnitude of the 
nuclear k^- enhancement may be obtained by compari- 
son of the p+Pb and p+C spectra. The spectra are 
compared to 7r° spectra measured in Pb+Pb collisions 
at y / SNN = 17.3 GeV with the same WA98 experimental 
setup [l6j |. 



In the WA98 experiment tt° yields were measured by 
detection of photons from the tt° — > 77 decay branch with 
a highly-segmented lead-glass calorimeter. This detector 
was located 21.5m downstream from the target and sub- 
tended the pseudorapidity range 2.3 < ?yi a b ^3.0. A 400 
GeV/c proton beam from the CERN SPS impinged on 
a beryllium production target to provide a mixed sec- 
ondary beam selected to have momentum of 160 GeV/c. 
The secondary beam consisted primarily of protons and 
pions with roughly equal content. Protons were iden- 
tified with use of two gas Cherenkov counters located 
upstream of the 1879 mg/cm 2 12 C (495 mg/cm 2 208 Pb) 
target. The WA98 minimum bias trigger condition re- 
quired a minimum amount of transverse energy (£+) in 



the region 3.5 



< 



?7lab 



< 



5.5, measured with a sampling 



calorimeter with electromagnetic and hadronic sections. 
The measured minimum bias cross section a m b for p+C 
(p+Pb) of 193 mb (1422 mb) corresponds to 86% (81%) 
of the total geometric cross section. The number of ana- 
lyzed minimum bias events was 1.2-10 6 (1.0T0 6 ) for p+C 
(p+Pb). A high-energy photon (HEP) trigger based on 
the sum energy signal of overlapping 4x4 groups of tow- 
ers in the lead-glass calorimeter was used to enhance the 
sample of high px events. The efficiency of the HEP trig- 
ger reached ~ 100% for photons with px ^ 0.8 GeV/c. 
An additional 1.5-10 6 (0.5-10 6 ) p+C (p+Pb) HEP events 
were analyzed which corresponded to 3.9 • 10 7 (8.2 • 10 6 ) 
sampled minimum bias events. 

Neutral pion yields were determined statistically by 
counting photon pairs with invariant mass in the tt° 
mass range after subtraction of the normalized back- 
ground from uncorrelated pairs. The shape of this back- 
ground was determined by mixing photons from different 
events. Only photon pairs with an energy asymmetry 
a = \Ei — E 2 \/ {Ei + E 2 ) < 0.7 were used in the analysis. 
A correction for geometrical acceptance and reconstruc- 
tion efficiency was applied to the raw ir° yields. The 
reconstruction efficiency takes into account the loss of 
7r s due to the photon identification and energy asym- 
metry cuts. It also modifies the ir° yield to account for 
the px shift that results from the finite energy resolution 
of the lead-glass calorimeter convoluted with the steeply 
falling 7T° px spectrum. Effects of overlapping show- 
ers in the calorimeter, which were important in central 
Pb+Pb collisions, are negligible in p+C and p+Pb colli- 
sions. The dominant systematic uncertainties are listed 
in Table |U The systematic uncertainties of the peak ex- 
traction, acceptance correction, and efficiency correction 
are approximately independent of px- The energy scale 
of the calorimeter was confirmed by comparison of the 
measured px-dependent ir° peak positions with GEANT 
simulations. The estimated uncertainty of 1.5% on the 
energy scale leads to an uncertainty on the 7r° yields that 
increases with px. 

The spectra of the invariant ir° yields in p+C and 
p+Pb collisions at -^/snn = 17.4 GeV are shown in Fig.[T] 
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TABLE I: Systematic uncertainties (in %) on the tv° yields in 
p+C and p+Pb collisions at three pt values. 
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FIG. 1: Invariant 7r° yields in minimum bias p+C and p+Pb 
collisions at ^/snn = 17.4 GeV. The error bars represent the 
quadratic sum of statistical and systematic uncertainties. The 
lines represent next-to- leading-order pQCD calculation of the 
7T° yield in y^s = 17.4 GeV p+p collisions for three different 
scales, scaled by (A^ ) = 1.7. The inset shows a comparison 
of the p+C spectrum to a parameterization of the yield in p+p 
collisions at ^/s — 17.3 GeV from [gjj. The dark gray box 
indicates the uncertainty of the parameterization estimated 
in [lj|] 5 the light gray box the uncertainty of (A^ ). 



a reliable reference at CERN SPS energies. The inset 
of Fig. [T] shows that the ir° yields per nucleon-nucleon 
collision measured in p+C are in good agreement with a 
parameterization of tt° spectra in p+p from Blattnig et 
al. [2(| that has been employed to study nuclear effects 
in Pb+Pb collisions [Hj]. 

Nuclear effects in ir° production can be quantified with 
the nuclear modification factor defined as 



(iV c P+ B ) dN^/d Pl 



IA+A 



{N A+A, 



"coll 



dA^o /dp] 



(1) 



Ip+B 



In the absence of nuclear effects -R^a ^ s expected to be 
unity for pt > 2 GeV/c where hard scattering is ex- 
pected to dominate particle production. (N co u} was de- 
termined with a Glauber Monte Carlo calculation [24| us- 
ing an inelastic nucleon-nucleon cross section of = 
(31.8 ± 2) mb (25|. The same Glauber calculation was 
used to extract (N co \\) values in Pb+Pb collisions. In 
the Glauber calculation the transverse energy Et was 
modelled by sampling a negative binomial distribution 
to determine the Et contribution of each participating 
nucleon 24| . The bias due to the trigger selection on the 
measured Et was taken into account. For minimum bias 
p+C and p+Pb collisions values of (N co \\) p+ c = 1.7±0.2 
and (Ac ii) p+ pb = 3.8 ± 0.4 were obtained. 

The (A^oii) values for Pb+Pb collisions were deter- 
mined by applying cuts to the simulated Et that cor- 
responded to the same fraction of cr, 1 ^^ 13 as the cuts 
applied to the measured Et. These (N co \\) values are 
listed in Table [TT] and agree within systematic uncertain- 
ties with those of Ref. [l6| that were determined with the 
VENUS 4.12 event generator in which cr?™ = 29.6 mb 
was used. With the large acceptance of the WA98 Et 
measurement, and good description of the Et distribu- 



tion, including fluctuations |14(, a centrality class corre 



sponding to the 1 % most central Pb+Pb collisions could 
be defined to access very large (N co u) values. 



The transition between the minimum bias and the HEP 
sample occurs at pt = 1.7 GeV/c. The p+C spectrum is 
compared to a next-to-leading-order (NLO) pQCD cal- 
culation for p+p at the same energy [21[ . The calculation 
was performed with the CTEQ6M [22j parton distribu- 
tion functions and the "Kniehl-Kramer-Potter" (KKP) 
set of fragmentation functions [23} with renormalization 
and factorization scales set to be the same at /i = Pt/2, 
Pt, or 2pT- For comparison with the p+C data the pQCD 
calculations have been scaled by the average number of 
nucleon-nucleon collisions in p+C. The pQCD calcula- 
tion shows a large uncertainty related to the arbitrary 
choice of scale. It has been suggested that the NLO per- 
turbative expansion is not sufficient at low energies and 
that threshold resummation corrections must be taken 
The large theoretical uncertainties 



21] 



into account 

demonstrate that pQCD calculations cannot be used as 



TABLE II: Results of the Glauber calculation for Pb+Pb col- 
lisions at y / SNN = 17.3 GeV (c™ = 31.8 mb). Centrality 
classes are given as a fraction (%) of o"^ +Pb ~ 6300 mb [l6| ]. 
The systematic uncertainty of (A co n) was parameterized as 
<4 coll = 0.1 2 + (0.4exp(-A coll /100)) 2 . 

class 123 4 5 6 7 86-8 

from 82.8 67.0 48.8 25.3 13.0 6.8 1.0 

to 100 82.8 67.0 48.8 25.3 13.0 6.8 1.0 13.0 

(A part } 8.2 23 54.2 123.2 218.2 289.1 347.8 382.6 322.5 

(Aeon) 6.3 22.1 67.1 202.9 433.1 627.0 803.7 912.0 727.8 

The p+Pb 7T° spectrum appears to be flatter than 
the p+C spectrum (see Fig. [T]). The ratio of the A^oii- 
normalized p+Pb and p+C spectra is shown in Fig. [2] 
At low pt (~ 1 GeV/c) the ratio is consistent with scal- 
ing with the number of participating nucleons (A^ part = 
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FIG. 2: The A^ co n-normalized ratio of 7r° yields in p+Pb and 
p+C collisions. The box at R' 



pPb 



0.8 represents the ex- 
pectation for the case of scaling with A/part rather than N co n- 
The box at unity indicates the N co n systematic uncertainty. 
The error bars represent the quadratic sum of the statisti- 
cal and remaining systematic uncertainties. The energy scale 
uncertainty was assumed to cancel. 
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FIG. 3: 7T° i?AA in Pb+Pb collisions at v /s^" = 17.3 GeV 
for three centrality classes using p+C or p+Pb spectra as 
a reference. The boxes around unity reflect the systematic 
uncertainties related to (N co n). 



N co \\ + 1 in p+A) while scaling with iV co n is observed at 
higher px (~ 2 GeV/c). In fact, forpx ^ 2 GeV/cthe ra- 
tio tends to be above unity, consistent with an expected 
stronger nuclear fcx-enhancement in p+Pb than in p+C. 

Neutral pion spectra from Pb+Pb collisions at ^/snn = 
17.3 GeV were published by WA98 in 0. Fig. 1 shows 
the 7T° nuclear modification factor -R^A as defined in Eq.Q] 
for three Pb+Pb collision centralities. For the 25.3 — 
48.8 % (of CTjnb +Pb ) most central collisions the data points 
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FIG. 4: Centrality dependence of the average (R'aa) f° r 2 < 
p T < 2.5 GeV/c and 2.5 < p T < 3.0 GeV/c. The bands 
around unity indicate the systematic uncertainty of (N co u). 
The error bars represent the quadratic sum of the statistical 
and systematic uncertainties. The jVp ar t values for the full 
points were shifted by +5 for better visibility. 



at highpx (pt ^ 2 GeV/c) suggest a stronger nuclear kr- 
cnhancement than in p+C. On the other hand, for the 
0—13% most central Pb+Pb collisions i?^A ^ s smaller 
than for the 25.3 — 48.8 % class but still consistent with 
unity at high px • For the — 1 % most central collisions 
the 7T° yield is significantly suppressed compared to either 
p+Pb or p+C as reference. An apparent suppression of 
the 7T° yield in very central Pb+Pb collisions compared 
to peripheral collisions was noted in Ref. (lq |. However 
the large uncertainty of using the peripheral distribution 
as a reference did not allow to draw a firm conclusion. 
The observed suppression is qualitatively consistent with 
expectations from jet-quenching. 

The centrality dependence of the average R' AA for 
2 < p T < 2.5 GeV/c and 2.5 < p T < 3.0 GeV/c is shown 
in Fig.[U The tt° yields in Pb+Pb in these px ranges are 
not suppressed for -/Vp art < 300 {(R'aa) ~ !)• F° r more 
central Pb+Pb collisions (R'aa) decreases with central- 
ity indicating significant suppression of the high px t° 
yield. The apparent lack of suppression, or enhancement 
even, for < 300 may be due to competing effects 

of suppression due to parton energy loss and nuclear k^- 
enhancement. 

In summary, 7r° spectra were measured in minimum 
bias p+C and p+Pb collisions at ^/snn = 17.4 GeV in 
the range 0.7 < px ^ 3.5 GeV/c. Based on these spectra 
the nuclear modification factors R'aa f° r Pb+Pb colli- 



sions at CERN SPS energies could be determined using 
a measured p+A reference. In very central Pb+Pb col- 
lisions (0 — 1 % of (Tfnb +Ph ) a significant suppression of 
high-px neutral pions was observed (R'aa ~ 0.5 ± 0.14) 
that is reminiscent of the high-px hadron suppression 
observed in Cu+Cu and Au+Au collisions at RHIC. The 
pion suppression reported here, together with the results 
at higher energies from RHIC, constrain jet-quenching 
models and may help clarify the behavior of the energy 
loss of partons in strongly interacting matter. 
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